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Abstract. Recent experimental breakthroughs in the treatment of dilute 
Bose gases have renewed interest in their quantum mechanical description, 
respectively in approximations to it. The ground state properties of dilute 
Bose gases confined in external potentials and interacting via repulsive short 
range forces are usually described by means of the Gross-Pitaevskii energy 
functional. In joint work with Elliott H. Lieb and Jakob Yngvason its sta- 
tus as an approximation for the quantum mechanical many-body ground state 
problem has recently been rigorously clarified. We present a summary of this 
work, for both the two- and three-dimensional case. 



1. Introduction 

The Gross-Pitaevskii (GP) functional was introduced in the early sixties as a 



phcnomenological description of the order parameter in superfluid He4 [G61, P61 
G63 1 . It has come into prominence again because of recent experiments on Bose 



Einstein condensation of dilute gases in magnetic traps. The paper [DGPS99| brings 
an up to date review of these developments. 

The present contribution is based on the joint work ]LSY00a| , |LSY00b| with 
Elliott H. Lieb and Jakob Yngvason (see also [LSYOOcQ. The starting point of our 
investigation is the Hamiltonian for N identical bosons moving in R 13 , D = 2 or 
3, that interact with each other via a radially symmetric pair-potential v{\xi — Xj\) 
and are confined by an external potential V(x): 



N 



(1) 



H = £{-A, + V(Si)} 



l<i<j<N 



The Hamiltonian acts on symmetric wave functions in <& N L 2 (M. D , dx). The pair 
interaction v is assumed to be nonnegative and of short range, more precisely, we 
demand it to have a finite scattering length. (For a definition of the scattering 
length in arbitrary dimension see | LY00 |.) The potential V that represents the 
trap is locally bounded and V(x) — > oo as \x\ — > oc. By shifting the energy scale 
we can assume that min^ V(x) = 0. 

Units are chosen so that h = 2m = 1, where m is the particle mass. A natural 
energy unit is given by the ground state energy fkv of the one particle Hamiltonian 
— (H 2 /2m)A + V. The corresponding length unit, y/ h/(jnui), measures the effective 
extension of the trap. 
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We are interested in the ground state energy E°- M = inf speciJ. Besides ./V it 
depends on the potentials V and v, but with V fixed and 

(2) v(r) = {a 1 /a) 2 vi(a 1 r/a), 

where v% has scattering length a\ and is regarded as fixed, E°- M is a function of N 
and a only. The corresponding ground state density is given by 

(3) p QM (x) = N [ \y a (x,x 2 ,...,x N )\ 2 dx 2 ...dx N , 



where \l/o is a ground state wave function of H. 

Note that v given in (||) has scattering length a. Here a is dimensionless and really 
stands for a^mui/h. Hence a scaling of v like (|J) is equivalent to scaling the external 
potential V at fixed v. In particlar the limit a — > with fixed V is equivalent to 
the limit u> — > with fixed v, if one introduces the scaling V(x) = wTq(w 1//2 x) for 
some fixed V\. 

Recent experiments on Bose-Einstein condensation are usually interpreted in 
terms of a function $ GP (:r) of x G K 13 , which minimizes the Gross-Pitaevskii energy 
functional 

(4) £ GP m= [ (|V$| 2 + ^|$| 2 +4vr.g|$| 4 )df 

under the subsidiary condition J |<I>| 2 = N. The corresponding energy is 

(5) E GP (N,g) = inf £ GP [$] = £ GP [$ GP ]. 

/ |$| 2 =AT 

The parameter g is different in dimensions 2 and 3. However, for any v alue of g > 
and N > it can be shown that a unique, strictly positive 3? exists [LSYOOa]. It 



depends on these parameters, of course, and when this is important we denote it 

-N,g- 



by $ GP 



The motivation of the term 47rg|$| 4 in the GP functional comes from the ground 
state energy density, £o(p), of a a dilute, thermodynamically infinite, homogeneous 
Bose gas of density p, interacting via a repulsive potential with s cattering length a 



The formulas for this quantity are older than the GP functional B47, HY57, S71|, 



at least for D = 3, but they have only very recently been derived rigorously for 



suitable interparticle potentials. See LY9£ | and LYOC ] . They are given by 



s (p) « 4:Trap z for D = 3, 

(6) e (p) w 47rp 2 |ln(a 2 j o)|- 1 for D = 2, 

where w means that the formulas are valid for dilute gases, where aP p <C 1. Hence 
the natural choice of the parameter g is 

(7) g = a for D = 3, 

(8) g = | ln(a 2 ^) | for£> = 2, 
where p is the mean GP density 

1 ^ 1 ^GP/ 



(9) p=-J |^ p (x)| 4 df. 

Note that <I> GP depends on g, so (0) together with (||) are non-linear equations 
for g. Alternatively, one could define g using the minimizer for g = 1 in the 
definition of p. Since p appears only under a logarithm, this would not effect our 
leading order calculations. For the same reason one could use the TF minimizer 
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(see below) instead of the GP minimizer to define g. Note also that unlike in the 
three-dimensional case g depends on N in the two-dimensional case. 

The idea is now that with this choice of g one should, for dilute gases, have that 

(10) E GP ^E°- M and p°- M (x) a |$ GP (x)| 2 = p GP (x). 
This is made precise in the following theorems. Note that by scaling 

(11) E GP (N,g)=NE GP (l,Ng) and ^ GP g (x) = N l /H GP Ng (x). 

Hence Ng is the natural parameter in GP theory. With this in mind we can state 
our first main result. 

Theorem 1.1 (The GP limit of the QM ground state energy and density). 

If N — > oo with Ng fixed, then 

, . , £Q M (iV,a) 

12) lim ' { = 1, 

V ; N^oo E GP {N,g) 

and 

(13) lim ^P QM (S) = |« 9 (^)| 2 

in the weak L\-sense. 

Note that by hypothesis of the theorem above it really applies to dilute gases, 
since for fixed Ng (which we refer to as the GP case) the mean density p is of order 
N and 

(14) a 3 p~N- 2 for D = 3, a 2 p ~ exp(-iV) for D = 2. 

Especially for D = 2 this is an unsatisfactory restriction, since a has to decrease 
exponentially with N. For a slower decrease Ng tends to infinity with TV, and 
the same holds for D = 3 if a does not decrease at least as iV -1 . In this case, the 
gradient term in the GP functional becomes negligible compared to the other terms 
and the so-called Thomas-Fermi (TF) functional 

(15) £ TF [p] = f (Vp + 47rgp 2 )dx 

arises. It is defined for nonnegative functions p on M. D . Its ground state energy 
E TF and density p TF are defined analogously to the GP case. Our second main 
result is that minimization of ( [l5| ) reproduces correctly the ground state energy 
and density of the many-body Hamiltonian in the limit when N — > oo, oP p — > 0, 
but Ng — > oo (which we refer to as the TF case) , provided the external potential is 
reasonably well behaved. We will assume that V is asymptotically equal to some 
function W that is homogeneous of some order s > and locally Holder continuous 



(see [LSYOObj for a precise definition). This condition can be relaxed, but it seems 



adequate for most practical applications and simplifies things considerably. 

Theorem 1.2 (The TF limit of the QM ground state energy and density). 
Assume that V satisfies the conditions stated above. If 7 = Ng — > 00 as N — > 00, 
but still a p — > 0, then 

(16) hm EQ 2 N ^ - 1, 

V 7 N^oo E TF (N,g) 
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and 

(17) lim ^^V M ( 7 1/(S+D) f) = p TF (f) 

N^QG iV 

in the weak Li-sense, where /5 TF is the minimizer of the TF functional under the 
condition J p = 1, g = 1, and with V replaced by W . 

Remark. The theorems are independent of the interaction potential v\ in (Q) . This 
means that in the limit we consider only the scattering length effects the ground 
state properties, and not the details of the potential. Note also that the particular 
limit we consider is not a mean field limit, since the interaction potential is very 
hard in this limit; in fact the term 47rg|$| 4 is mostly kinetic energy. 



In the following, we will present a brief sketch of the proof of Theorems 1.1 and 



1.2 . We will derive appropriate upper and lower bounds on the ground state energy 
fiQM _ rp-^g conver g ence f the densities follows from the convergence of the energies 
in the us ual w ay by va riation with respect to the external potential. We refer to 



|LSY00a| and flLSY00b| for details. 



2. Upper bound to the QM energy 

To derive an upper bound on E^ M we use a generalization of a trial wave function 
of Dyson [ D57 ], who used this function to give an upper bound on the ground state 
energy of the homogeneous hard core Bose gas. It is of the form 

N 

(18) y(x!,...,x N ) = Y[<S> GP (x l )F(x 1 ,...,x N ), 

i=l 

where F is constructed in the following way: 

N 

(19) F(x u . . .,x N ) = JI /(ti (a?i, ■ ■ ■ ,a?i)), 

i=l 

where = min{|a?i — Xj\, 1 < j < i — 1} is the distance of Xi to its nearest neighbor 
among the points aft, ... , a?j_i, and / is a function of r > 0. We choose it to be 



(20) /(r) = 



fo{r) /f {b) for r<b 
1 for r > 6, 



where /o is the solution of the zero-energy scattering equation (see [LY00Q and b is 
some cut-off parameter of order b ~ p~ x l D . The function is not totally sym- 
metric, but for an upper bound it is nevertheless an acceptable test wave function 
since the bosonic ground state energy is equal to the absolute ground state energy. 
The result of a somewhat lengthy computation is the upper bound 



(2i, ^ J v,«)< e - (w ,p)x{;+°(|'Q. )| _ p 



1 + Oiap 1 / 3 ) for L> = 3 

for D = 2. 



with the power p equal to 1 in the GP case and s/(s + 2) in the TF case (where V 
is asymptotically homogeneous of order s). 
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3. Lower bound to the QM energy 
To obtain a lower bound for the QM energy the strategy is to divide space into 



boxes and use the estimate on the homogeneous gas, given in [ LY9S ] and [ LY00 |, in 
each box with Neumann boundary conditions. One then minimizes over all possible 
divisions of the particles among the different boxes. This gives a lower bound to 
the energy because discontinuous wave functions for the quadratic form defined by 
the Hamiltonian are now allowed. We can neglect interactions among particles in 
different boxes because v > 0. Finally, one lets the box size tend to zero. However, 
it is not possible to simply approximate V by a constant potential in each box. To 
see this consider the case of noninteracting particles, i.e., v = and hence a = 0. 
Here E°- M = Nfkv, but a 'naive' box method gives only as lower bound, since it 
clearly pays to put all the particles with a constant wave function in the box with 
the lowest value of V. 

For this reason we start by separating out the GP wave function in each variable 
and write a general wave function ^ as 

N 

(22) ^f(S 1 ,...,x N ) = Y[^ GP (St)F{2 u ...,x N ). 

i=l 

This defines F for a given because $ GP is everywhere strictly positive, being the 
ground state of the operator — A + U+87r<7|$ GP | 2 . We now compute the expectation 
value of H in the state \P. Using partial integration and the variational equation for 
<£> GP , we see that for computing the ground state energy of H we have to minimize 
the normalized quadratic form 

(23) 
Q{F) 



-l 

— ^ 



N lUtl P GP (^) IV^I 2 + E vQZi - x 3 \)\F\ 2 - Sngp^ix^Ff 



^ /n£=i/> GP (**)m 2 



2=1 

Compared to the expression for the energy involving itself we have thus obtained 
the replacements 

N N 

(24) V(x) -> -8irgp GP (x) and JJdSi -> JJ p GP {x^dx, 

i=l i=l 

(recall that p (x) = |<f >GP (a f )| 2 ). We now use the box method on this problem. 
More precisely, labeling the boxes by an index a, we have 

(25) miQ(F) > inf VinfQ Q (F Q ), 

a 

where Q a is defined by the same formula as Q but with the integrations limited 
to the box a, F a is a wave function with particle number n a , and the infimum is 
taken over all distributions of the particles with ^]ii a — N. Approximating p GP 
by a constant in each box, we can use the bound on the homogeneous case ([ LY98| 



and [ LYOOj ] ) in each box. To control the error terms, we need a lower bound on 



the ratio p a . m i n / ' p a . max in each box, which stems from the measures in (|23|), where 
Pa, max and jQ a ,mm> respectively, denote the maximal and minimal values of p in 
box a. 
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The problem is that for any hxed size of boxes pa,min/Pa,max tends rapidly to 
zero for boxes far from the origin. This problem can be solved by enclosing the 
whole system in a big box A# of side length R, with Neumann conditions on the 
boundary. Replacing $ GP by , which is the minimizer of £ GP restricted to A# 
(and satisfies Neumann conditions), and restricting the integrations to we can let 
the side lengths of the small boxes tend to zero and be sure that p a ,min/7>a,max - * 1 
uniformly for all the boxes a. However, we must control the error made by enclosing 
the system in the big box. Let E^ M (N, a) denote the quantum mechanical ground 
state energy in the box A#. The essential step is 

Lemma 3.1. There is an Rq < oo, depending only on Ng, such that 

(26) E Q - M (N,a)>E% M (N,a) 

for all R > i?o o,nd all N , a with Ng fixed. 

This lemma follows from V{x) — > oo for \x\ — > oo, together with an estimate for 
the chemical potential 

(27) E% M (N + 1, a) - E% M (N, a) < e(Ng)(l + o(l)) 

where e(Ng) depends only on Ng and is independent of R. The proof of ( p7j ) is 
similar to the proof of the upper bound (|2l|). 

Now p a ,min/Pa,max is bounded below uniformly in each small box contained in 
Afl. We first let the size of the boxes tend to zero as N — > oo, and finally take the 
limit R — > oo. This proves the desired lower bound (for fixed Ng, i.e. the GP case). 
If Ng — > oo as N — > oo the method above does not work since the Rq of Lemma 



3.1 tends to infinity with Ng. However, using the explicit form of the TF minimizer, 



namely 



1 



(28) P T N F g (x) = —[p TF -V(x)} + , 

where [t]+ = max{i, 0} and p TF is chosen so that the normalization condition 
J pjf g = N holds, we can use V(x) > p TF — 8irgp TF (x) to get an replacement 
as in (|24|) without changing the measure. Moreover, p TF has compact support, so, 
applying again the box method described above, the boxes far out do not contribute 
to the energy. However, p TF (which depends only on the combination Ng) tends 
to infinity as Ng — » oo. Therefore we need to control the asymptotic behaviour 
of the external potential, and this leads to the restrictions on V described in the 



paragraph preceding Theorem 1.2 



After controlling all error terms properly, we arrive at the result 

(29) limM E rr J N,a } > 1 

V ; w^oo E TF {N,g) ~ 

in the limit N — > oo, aPp — > and Ng — > oo. Together with the upper bound 
© and the fact that E GP (N , g) / E TF {N , g) = E GP {1, Ng)/E TF {1, Ng) 1 as 
Ng — » oo this proves Theorem |l.2| . 
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